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INTRODUCTION 


accurate calculation of aeroelaatic characteristics re<,uir.d for the analysis and design of high- 
perfornance aircraft requires accurate and efficient evaluation of stead, «.d unsteady aerodynanic loads 
cn aircraft having arbitrary shapes and notions, including structural de'or«tions. This presentation 
uill address the aircraft geonetry re,uirenents for unsteady aerodyn»nic conput.tions and .ill enphasite 
differences between requirements for steady and unsteady flow. 
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COMPARISON OF STEADY AND UNSTEADY FLOW REQUIREMENTS 

(Figure 1) 

Requlretrents for aeroelastic analysis and design are in several respects more complicated and more 
severe than for the more convent l->nal steady-state aeroJynanics. For example: (1) The aeroelasticlan 

deals with flexible structures so that even in steady-state conditions, the aerodvnamic load is a function 
of structural deformation, and vice versa. (2) The unsteacy aerodynamic formulations required in dynamic 
aeroelast J city involve complex quantities (e.g.. normalwas i velocities, aerodynamic influence functions, 
and pressure) that manifest time- or frequency-dependent aLtenuations and phase shifts relative to steady* 
state. (3) In dynamic aeroelastic! ty— flutter . for example— the aeroelasticlan must evaluate pressure 
distributions for ^-Ibratlon mode shapes that are much more vlggly than a typical steady-state mean-camber 
surface. The corresponding pressure distributions will also be more wiggly than those for steady state 
so that coiLputatlonal convergence requirements are usually more severe than for steady state. (4) Flutter 
analyses, as well as iterative structural resizing, require evaluation of pressure distributions for a 
multiplicity of mode shapes, frequencies, aircraft loading conditions, etc. Consequently, computational 
efficiency is vital, and It Is essential to minimize the amount of recomputation required when mode shapes 
and/or frequencies are changed. 

With these thoughts In mind, we shall discuss geometry requirements within the framework of the 
SOUSSA aerodynamic formulation because It is the most general potential-flow program that we now have 
under development (with regard to aircraft geometry, motion and deformations, and speed ranges) and 
because present and future SOUSSA geometry requirements are as stringent as those for any aerodynamic 
program that we now anticipate. Geometry required is considered to be composed of three parts: (1) 

shape of vehicle. (2) orientation of vehicle. (3) deformat ion(s) . Orientation involves little more 
than a rotation of coordinate axes and consequently will not be emphasized here. Deformations can be 
finite but are more usually taken to be infinitesimal and approximated by a linear combination of the 
natural undamped vibration modes of the aircraft. As many as two dozen modes or more may be required 
to converge the aeroelastic solution. A related geometrical requirement is determination of wake 
shape which is not known a priori although it may be assumed to be flat for many applications* 



STEADY, OSCILLATORY AND UNSTEADY SUBSONIC AND SUPERSONIC AERODYNAMICS 

(SOUSSA) 


Objective; An accurate, general, unified* method for calculating steady and unsteady loads on 
complete aircraft with arbitrary shape and motion in subsonic or supersonic flow, 
with emphasis on application in computer-aided structural design 

Approach: Green's theorem is used to formulate exact Integral equation for potential'. 

where perturbation velocity potential 

G = Green's function 
F = nonlinear terms 
Sfx, y, z. t) = 0 defines body surfa ce 

Exact boundary condition on body is 


Figure 2 



STEADY, OSCILLATORY AIu) UNSTEADY SUBSONIC AND SUPERSONIC A£hODY:;A.MTCS 

(SOUSSA) 

(Figure 2) 

To set the stage, a brief review of SOUSSA formulation is in order. Application of Green’s theorem 
leads to an integral equation for the perturbation velocity potential ^ at any point P in the flew 
or on the flow boundary (i.e., on S). Note that the second integral contains only ^...near terms which 
are integrated over the boundary surface S, whereas the first integral contains nonlinear terms F, 
involving products of derivatives of ifi, vdiich must be integrated over the fluid volume. 

The boundary condition clearly shows the effect of time variation of S. If the variation is 

3S 

harmonic, for example, the ^ term becomes iuS so that the normal wash at S = 0 becomes complex. 

The imaginary part, however, involves only surface ordinates (including displacements and deformations), 
whereas the real (steady-<?tat») part involves derivatives of surface ordinates. Ilms, introduction 
of unsteadiness doos not impose more stringent requirements on surface definition as far as quantities 
required are concerned. It may, however, reqxiire greater accuracy and greater amounts of geometrical 
information (e.g. , for more points on the surface) in order to define at’equately the wiggly modes of 
deformation referred to previously. 

The influence of nonlinear terms F is being studied in the develoisaent of SOUSSA aerodynamics for 
the transonic range. However, these terms are not included in the present computer program. 
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SOLUTION OF INTEGRAL EQUATION 


cc 


Solution by spatial discretization with arbitrary nonplanar quadrilateral 
surface panels and time solution by Laplace transform results in 



where 


^ = Laplace transform of perlurbalioii velocity potential 
^ = Laplace transform of noniialwjsh 







s = Laplace transform variable 

B ju, C jh, Djh, F ;p, G m = integrals over surface panels, independent 

of normalwash and s 

0jl^, rijh = lag functions 


Surface pressures are obtained from Bernoulli's equation. 


Figure 3 



SOLUTION OF INTEGRAL EQUATION 


(Figure 3) 

AssuBning that the effect of the motion of the surface is negligible except in the boundary 
condition (infinitesimal unsteady deformation) and neglecting the nonlinear terms » then surface 
paneling and Laplace transform ikolution yield a matrix equation relating the unknown potential $ 
on the vehicle surface to the normal wash ip* Elements of the coefficient matrices are independent 
of normalwash (and hence deformation) and are simple functions of the Laplace variable s. For a given 
paneling arrangement, they depend only on Mach number. 

Use of arbitrary nonplanar quadrilateral panels permits matching nod<i of the aerodynamic panels to 
the nodes of a structural finite-element model, if desired, in order to use the nodal coordinates and 
calculated displacements directly without requirement for interpolation. In general, however, solution 
for the velocity-potential matrix requires the following geometrical input: (1) Coordinates of panel 

nodes usually obtained by interpulat ion (lofting) from aircraft shape information. (2) Time-dependent 
normalwash at control points which usually do not coincide with panel nodes. Normalwash involves 
coordinates and slopes obtained from aircraft shape plus orientation plus deformation. 

Note that increasing the number of deformation modes used involves only adding columns to the and 
4) matrices, and that updating the entire set of deformation modes, as in a structural design application. 
Involves only replacing the matrix. The Y and Z matrices are unaffected in either case. 

Surface pressures are obtained from Bernoulli's equation. Generalized aerodynamic forces, including 
aerodynamic coefficients and stability derivatives, are determined from weighted integrals of the pressure 
which require values of surface displacement (due to rigid-body rotation and/or i&odal deformation) at a 
set of integration points which may not coincide with the panel nodes nor control points. 

The geometrical information required by SOUSSA can, of course, be generated with any suitable 
geometry preprocessor as long as the results are cast in required SOUSSA input format. It is evident, 
however, that automatic paneling capability is essential to the efficient processing of complicated 
shapes and deformations that may require many hundreds cf panels. Such capability should include not 
only automatic calculation of the coordinates of nodes, control points, and integration points, but 
also automatic identification numbering for these points as well as for the panels and systematic 
identification of vdiich nodes go with which panels. 
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GENERAL POTENTIAL-FLOW AERODYNAMICS 

(SOUSSA) 

GENERAL FINITE -ELEMENT METHOD: 

©ARBITRARY COMPLETE A/C CONFIGURATION 

® STEADY AND GENERAL UNSTEADY MOTION 

® SUBSONIC AND SUPERSONIC 

o COMPUTATIONAL EFFICIENCY 

CURRENT DEVELOPMENTS: 

® NONLINEAR EFFECTS (TRANSONIC FLOW 
WAKE DEFORMATION) 

o IMPROVED FINITE ELEMENTS (HIGHER ORDER 
SPECIAL PURPOSE) 

0 ROTATIONAL FLOW (TURBULENCE. VISCOSITY) 


Figure 4 


^EliERAL POTEKTIAI/-FLOW AERODY!?yCCS 


(SOUSSA) 

(Figure U) 

The upper part of this figure lists some of the fea ores of the SOUSSA aerodynamic formulation 
The lower part indicates son* expanded capabilities and inprcveaents that are under development and 
that will influence surface geometry requirements. These are discussed in the following figures. 



SURFACE OF INTEGRA! iON S 



(f continuous 
Rank ine-Hugon iot 
Variable location 


Fixed location 


p continuous 


Figure 5 



SURFACE OF lOTEGRATION S 


(Figure 5) 


In the integral equation for the velocity potential (previously shown) the surface integration extends 
over the surface of the aircraft plus its wake, and the no- penetration boundary condition DS/Dt * O 
applies over both. Moreover, the pressure -'ust be continuous across the wake although the potential is 
discontinuous. Ihe forward edge of the wake of a lifting surface is considered to be located at the 
lifting-surface trailing edge, but the position of the rest of the wake is not known a priori, is varlsble, 
time- dependent, and must be determined in the calculation. This variability requires relocation and 
reorientation of the wake and its panels during the calculation, perhaps many times if the calculation 
is iterative. 

This figure also shows a shockwave which is isolated from the flow field by a portion of the surface 
S. Over this portion of S the no-penetratlon boundary condition must be replaced by Rankine-Hugoniot 
conditions or other shock conditions which quantify siu>ck-induced discontinuities in derivatives of 
the potential although the potential itself is continuous across the shock. These discontinuities make 
it desirable to have panel edges lie along the foot of the shock. But shock location, shape, extent, 
strength, and velocity relative to the vehicle surface are time dependent. Moreover, motion of finite 
amplitude— even small amplitude— that is needed to investigate limit— cycle aeroelastic response can 
lead to large- amplitude shock motion and even discontinuous shock location. Consequently, requiring 
panel edges to coincide with the foot of the shock can require extensive repanellng in the vicinity 
of the shock during calculations for unsteady motion. In contrast, nonlinear calculations for shock- 
free transonic flow require no repaneling and impose no special requirements for surface geometry. 
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WAKE-FUSELAGE INTERSECTION 



Figure 6 


WAKE-FUSELAGE INTERSECTION 


(figure 


Shocks and wakes from lifting surfaces impinge upon fuselages or other portions of the 
vehicle along lines that are time dependent. Because of :he discontinuities In potential or 
its derivatives at these implngeinent lines * it is cfcsirable that panel edges coincide with 
them. Hence, time dependent repaneling in these vicinities is also indicated. 

For simplicity in its development, ^he present S^USSA program contains zeroth-order 
(constant-potential) aerodynamic elements. However, it has been Intended from the beginning 
that the program would employ higher-order elements in order to reduce the number ot elements 
required to converge the solution. Such elements have been developed and will soon be 
incorporated into the program. In addition, special-purpose elements are being develooed 
for paneling In regions where correct variation of potential is theoretically known. These 
elements have built-in shape functions to produce the correct variation of potential, for 
example, adjacent to normalwash discontinuities sucl» as control-surface hinge lines and side 
edges, or correct variation of potential derivatives ns at subsonic trailing edges. In addition, 
flow-through elements are required tc model engine thrust in nacelles and tc p.anel shockwaves. 
Such elements impose no new requirements for surface geosietry information. Note that no special 
panels arc required adjacent to shock or wake impingement on the body surface. Use of higher- 
order and special-purpose elements should reduce computer tisic and storage requirements but 
probably will do little to reduce the amount of geometrical input information required. 

Although fewer elements are used, more information is required per eleirent. Detailed accuracy 
of information out requl^w-s detailed accuracy of information In, irrespective of the level of 
sophistication. 

Finally, incorporating the effects of viscosity and rotational flc4ir will impose a 
requirement for relatively high accuracy of computed pressure gradients and hence will require 
higher-order elements (at least third order) than would be required for most potential-flow 
problems. Alternatively, it is possible that required accuracy and order of continuity may be 
attained from solutions usi/.g lover-order elements followed by spline (or other) interpolation of 
the calculated potential. 

Specific treatment of aircraft geometry in the SOUSSA program itself will next be described. 
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POTENTIAL-NORMALWASH RELATIONSHIP 
(Fi»»ure 7) 

The rest of this presentation Is devoted to the specific geometry 
requirements for the program SOUSSA P (Steady, Oscillatory and Unsteady, 
Subsonic and Supersonic Aerodynamics; P^roductlon Version). As presented 
above, the Green's function method yields an integral equation over the 
surface of the aircraft and its wake (with differential-delay dependence 
on time). Dividing the surfaces in quadrilateral elements and assuming 
the potential, the normalwash and the potential discontinuity to be 
constant within each element, one obtains Eq. (l). 

The coefficients B,. , C,. , etc., are evaluated analytically, with 

Jh jh 

the original siu*face approximated by a hyperboloidal paraboloid 

(hyperboloidal elem*'nt). Low-order numerical quadrature is used for 
distant elements. 

In order to complete the formulation three additional relationships 
are required; 

1. Boundary conditions, relating normalwash to the generalized 
coordinates q^ (Eq. 2) 

2. Bernoulli's theorem relating pressure coefficient to 

potential (Eq. 3) 

3. Definition of generalized forces, e , as functionals of the 

m 

pressure coefficient Cp (Eq. 1*) 

FinaQly combining Eqs. 1 to U one obtains the matrix E relating 
the generalized forces, e^, to the generalized coordinates, q^. 





















SOUSSA P 1.1 


(Figure 8) 
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CHECKPOINT #1 


(Figure 9) 

This figure preecnts the contents of Checkpoint #1 (input to module 
BODYO), which consists of inforrantlon describing the geometry of the 
aircraft body. It Is user-oriented in that the quantities required are 
compatible with the output of state-of-the-art geometry preprocessors. 
Also, if the aircraft is synuketrlc with respect to the x-z plane, then 
only the right half need be supplied. The same is true for the x-y plane. 
Regarding the individual components of Checkpoint #1: 
o The Cai-tesian coordinates of the nodes are assumed to be those of 
the aircraft already oriented as desired by the user, 
o Referring to the example depicted in the figure, element number 1, 
corner 1 yields » ode number 2. 

o The body-symmetry code numbers reflect whether symmetry is 
considered with respect to the x-z and/or y-z planes, 
o The element code numbers provide information such as whether 
or not a wake emanates from an edge of an element, or if an edge 
coincides with a hinge line, etc. 
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CHECKPOINTS #2, 3 AND »4 


(Figure 10) 

This figure presents the contents of Checkpoints 02 (input to Module 
WAKEG), 03 (input to Module EUMOD) and 0k (input to Module ElMOD). 

Checkpoint 0Z consists of infomation describing the geometry of the 
vake. Dy describing the wake as a collection of strips, many different 
forma of input can easily be made compatible. If a wake strip is symmetric 
with respect to the x-z plane, then only the right half need be supplied 
(same for the x-y plane). 

Also, 

o Desired orientation of the wake with respect to the aircraft is 
assumed to already be satisfied. 

o The matrix which relates each wake-strip number with the 
corresponding four trailing-edge element numbers is used in 
evaluating the trailing edge values of the potential and for 
determining the values of the pressure discontinuity at the 
centroid of the trailing-edge elements. 

Checkpoint 03 corresponds to the general! zed- forces deformation modes, 
and Checkpoint 0k corresponds to the boundary-condition deformation modes. 
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Figure 11 



CHECKlVlNrJ #5 AND 6 


(Figure 11) 

This figure presents the contents of Checkpoints (input to module COEFB) 
and #6 (input to modules COEFB and COKi’’W). These check{>oint8 are at a 
lower level tlian checkpoints #l-k if one views the SOUSSA P flow diagram 
as a top-down representation. The implications of this are that these 
checkpoints arc not as "user-oriented'' as higlier-level checkpoints, since 
program execution has progressed to this point. This is evidenced by the 
fact tliat for Checkpoint #5, the some quantities as Checkpoint 01 are 
required except that symmetry conditions (and their advantages in preparing 
geometrical input) are not considered. Furthermore, geometrical quantities 
such as the base vectors and normals of surface elements are not as readily 
available from geometry preprocessors as the information contained in 
Checkpoint 01 , These considerations must be accounted for by those 
users desiring co begin execution of SOUSSA P at this level. 

For version 1.1 of SOUSSA P, the location of the control points must 
be specified as the geometrical centroids of the body elements. For 
future versions (first-order finite element formulation), the location 
of the nodes will be the necessary input. 
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CHECKPonrr #7 


(Pif^ure 12) 

This figure presents the contents of Checkpoint #7 (input to nodule 
COEFW). This is also not a top-level checkpoint, hence, its contents siny 
not be as **u»er-oriented" as, say. Checkpoint 02 . For Instance, at this 
level : 

o The coordinates of the wake elements (as opposed to the wake 
strips) are required, and no symmetry conditions nay be taken 
advantage of (i.e. , all the elements must be input). 

o The matrix used in correction for the trai ling-edge potential 
values, for example, must be given for the elements comprising 
the complete woke. 

o Most geometry preprocessors would not provide the matrix of the 
coefficients of influence of the trailing-edgc elements that 
determines the value of the potential discontinuity for each woke 
element. Note for SOUSSA P 1.1 these coefficients are simply 1 
and -1, but for later versions, splines could be used to determine 


these coefficients. 
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Figure 13 
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(Figure 13) 

Work now under consideration Includes higher-order finite-element 
formulations. In particular, work on first-order finite element is now 
underway. Splines are also being considered for various nodules. 

Also under consideration is the extension to transonic unsteady flow. 
This requires the addition of two terms, the first is due to the r.oii- 
linearlties in the differential equation, the second to the presence 
of shocks. The equation is solved step-by-step in the time domain. The 
location of the shock and its intensity are obtained at each time stop. 


